Exponential Integral

In mathematics, the exponential integral Ei is a special function on the complex plane.
It is defined as one particular definite integral of the ratio between an exponential function and its argument.

Definitions

For real non-zero values of ''x'', the exponential integral Ei(''x'') is defined as

:$\operatorname(x) = -\int_^\infty \fract\,dt = \int_^x \fract\,dt.$

The Risch algorithm shows that Ei is not an elementary function. The definition above can be used for positive values of ''x'', but the integral has to be understood in terms of the Cauchy principal value due to the singularity of the integrand at zero.

For complex values of the argument, the definition becomes ambiguous due to branch points at 0 and Instead of Ei, the following notation is used,

:$E_1(z) = \int_z^\infty \frac\, dt,\qquad, (z), <\pi$

For positive values of ''x'', we have

In general, a branch cut is taken on the negative real axis and ''E''₁ can be defined by analytic continuation elsewhere on the complex plane.

For positive values of the real part of $z$, this can be written
:$E_1(z) = \int_1^\infty \frac\, dt = \int_0^1 \frac\, du ,\qquad \Re(z) \ge 0.$

The behaviour of ''E''₁ near the branch cut can be seen by the following relation:

:$\lim_ E_1(-x \pm i\delta) = -\operatorname(x) \mp i\pi,\qquad x>0.$

Properties

Several properties of the exponential integral below, in certain cases, allow one to avoid its explicit evaluation through the definition above.

Convergent series

For real or complex arguments off the negative real axis, $E_1(z)$ can be expressed as

:$E_1(z) = -\gamma - \ln z - \sum_^ \frac \qquad (\left, \operatorname(z) \ < \pi)$

where $\gamma$ is the Euler–Mascheroni constant. The sum converges for all complex $z$, and we take the usual value of the complex logarithm having a branch cut along the negative real axis.

This formula can be used to compute $E_1(x)$ with floating point operations for real $x$ between 0 and 2.5. For $x > 2.5$, the result is inaccurate due to cancellation.

A faster converging series was found by Ramanujan:

:$(x) = \gamma + \ln x + \exp \sum_^\infty \frac \sum_^ \frac$

These alternating series can also be used to give good asymptotic bounds for small x, e.g. :

:$1-\frac\le (x) - \gamma - \ln x \le 1-\frac+\frac$
for $x\ge 0$.

Asymptotic (divergent) series

Unfortunately, the convergence of the series above is slow for arguments of larger modulus. For example, more than 40 terms are required to get an answer correct to three significant figures for $E_1(10)$. However, for positive values of x, there is a divergent series approximation that can be obtained by integrating $x e^x E_1(x)$ by parts:
: $E_1(x)=\frac x \left(\sum_^ \frac +O(N!x^) \right)$
The relative error of the approximation above is plotted on the figure to the right for various values of $N$, the number of terms in the truncated sum ($N=1$ in red, $N=5$ in pink).

Exponential and logarithmic behavior: bracketing

From the two series suggested in previous subsections, it follows that $E_1$ behaves like a negative exponential for large values of the argument and like a logarithm for small values. For positive real values of the argument, $E_1$ can be bracketed by elementary functions as follows:
:$\frac 1 2 e^\,\ln\!\left( 1+\frac 2 x \right) < E_1(x) < e^\,\ln\!\left( 1+\frac 1 x \right) \qquad x>0$

The left-hand side of this inequality is shown in the graph to the left in blue; the central part $E_1(x)$ is shown in black and the right-hand side is shown in red.

Definition by Ein

Both $\operatorname$ and $E_1$ can be written more simply using the entire function $\operatorname$ defined as
:$\operatorname(z) = \int_0^z (1-e^)\frac = \sum_^\infty \frac$
(note that this is just the alternating series in the above definition of $\mathrm_1$). Then we have
:$E_1(z) \,=\, -\gamma-\ln z + (z) \qquad \left, \operatorname(z) \ < \pi$
:$\operatorname(x) \,=\, \gamma+\ln - \operatorname(-x) \qquad x \neq 0$

Relation with other functions

Kummer's equation
:$z\frac + (b-z)\frac - aw = 0$
is usually solved by the confluent hypergeometric functions $M(a,b,z)$ and $U(a,b,z).$ But when $a=0$ and $b=1,$ that is,
:$z\frac + (1-z)\frac = 0$
we have
:$M(0,1,z)=U(0,1,z)=1$
for all ''z''. A second solution is then given by E₁(−''z''). In fact,
:$E_1(-z)=-\gamma-i\pi+\frac,\qquad 0<(z)<2\pi$
with the derivative evaluated at $a=0.$ Another connexion with the confluent hypergeometric functions is that ''E₁'' is an exponential times the function ''U''(1,1,''z''):
:$E_1(z)=e^U(1,1,z)$

The exponential integral is closely related to the logarithmic integral function li(''x'') by the formula
:$\operatorname(e^x) = \operatorname(x)$
for non-zero real values of $x$.

Generalization

The exponential integral may also be generalized to

:$E_n(x) = \int_1^\infty \frac\, dt,$

which can be written as a special case of the upper incomplete gamma function:

: $E_n(x) =x^\Gamma(1-n,x).$

The generalized form is sometimes called the Misra function $\varphi_m(x)$, defined as

:$\varphi_m(x)=E_(x).$

Many properties of this generalized form can be found in th
NIST Digital Library of Mathematical Functions.

Including a logarithm defines the generalized integro-exponential function
:$E_s^j(z)= \frac\int_1^\infty \left(\log t\right)^j \frac\,dt.$

The indefinite integral:
:$\operatorname(a \cdot b) = \iint e^ \, da \, db$
is similar in form to the ordinary generating function for $d(n)$, the number of divisors of $n$:
:$\sum\limits_^ d(n)x^ = \sum\limits_^ \sum\limits_^ x^$

Derivatives

The derivatives of the generalised functions $E_n$ can be calculated by means of the formula
:$E_n '(z) = - E_(z) \qquad (n=1,2,3,\ldots)$
Note that the function $E_0$ is easy to evaluate (making this recursion useful), since it is just $e^/z$.

Exponential integral of imaginary argument

If $z$ is imaginary, it has a nonnegative real part, so we can use the formula
:$E_1(z) = \int_1^\infty \frac t \, dt$
to get a relation with the trigonometric integrals $\operatorname$ and $\operatorname$:
:
E_1(ix) = i\left -\tfrac\pi + \operatorname(x)\right- \operatorname(x)
\qquad (x > 0)

The real and imaginary parts of $\mathrm_1(ix)$ are plotted in the figure to the right with black and red curves.

Approximations

There have been a number of approximations for the exponential integral function. These include:

* The Swamee and Ohija approximation $E_1(x) = \left (A^+B \right )^,$ where \begin
A &= \ln\left left (\frac+0.65\right)(1+x)\right\\
B &= x^4e^(2+x)^
\end
* The Allen and Hastings approximation $E_1(x) = \begin - \ln x +\textbf^T\textbf_5,&x\leq1 \\ \frac x \frac,&x\geq1 \end$ where \begin
\textbf & \triangleq 0.57722, 0.99999, -0.24991, 0.05519, -0.00976, 0.00108T \\
\textbf & \triangleq .26777,8.63476, 18.05902, 8.57333T \\
\textbf & \triangleq .95850, 21.09965, 25.63296, 9.57332T \\
\textbf_k &\triangleq ^0,x^1,\dots, x^kT
\end
* The continued fraction expansion $E_1(x) = \cfrac.$
* The approximation of Barry ''et al.'' E_1(x) = \frac\ln\left +\frac G x -\frac\right where: $\begin h &= \frac+\frac \\ q &=\fracx^ \\ h_ &= \frac \\ b &=\sqrt \\ G &= e^ \end$ with $\gamma$ being the Euler–Mascheroni constant.

Applications

* Time-dependent heat transfer
* Nonequilibrium groundwater flow in the Theis solution (called a ''well function'')
* Radiative transfer in stellar and planetary atmospheres
* Radial diffusivity equation for transient or unsteady state flow with line sources and sinks
* Solutions to the neutron transport equation in simplified 1-D geometries

See also

* Goodwin–Staton integral
* Bickley–Naylor functions

Notes

References

*
Chapter 5

*
*
*
*
*
*
*
*
*
*
*
*
*

External links

*

NIST documentation on the Generalized Exponential Integral
*
*
*

Exponential, Logarithmic, Sine, and Cosine Integrals
in DLMF.

{{DEFAULTSORT:Exponential Integral
Exponentials
Special functions
Special hypergeometric functions
Integrals